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Abstract 

The short-distance coefficients for the color-octet 3 Si term in the fragmentation function for a 
gluon to split into polarized heavy quarkonium states are re-calculated to order a^. The light-cone 
gauge remarkably simplifies the calculation by eliminating many Feynman diagrams at the expense 
of introducing spurious poles in loop integrals. We do not use any conventional prescriptions for 
spurious pole. Instead, we only use gauge invariance with the aid of Collins-Soper definition of the 
fragmentation function. Our result agrees with a previous calculation of Braaten and Lee in the 
Feynman gauge, but disagrees with another previous calculation. 
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I. INTRODUCTION 



Heavy quarkonium state is one of the simplest hadron with which we can probe both 
perturbative and nonperturbative nature of quantum chromodynamics. Among various 
quarkonia, S-wave spin-triplet states are especially interesting because of their clean leptonic 
decay modes. Based on factorization theorems for inclusive single-hadron production Q, 
Q], one can deduce that dominant production mechanism for heavy quarkonia with large 
transverse momentum pt is fragmentation |3[ , the production of a parton which subsequently 
decays into the quarkonium state and other partons. In the region, where fragmentation 
dominates, theoretical calculation becomes simplified and more reliable. Measurements of 
large-pr S'-wave spin-triplet charmonia production cross section at the Fermilab Tevatron 
has led to a remarkable progress in heavy quarkonium physics based on nonrelativistic 
QCD (NRQCD) In high-energy pp collisions, gluon production rate is dominant and 
inclusive production of a heavy quark pair QQ via subsequent decay of this almost on-shell 
gluon is enhanced by the gluon propagator [5|. Unexpectedly large measured production rate 
of direct J / ip and ip' at large pt at the Tevatron (f| was explained by the gluon fragmenation 
into a color-octet QQ pair followed by a nonperturbative NRQCD transition into the spin- 
triplet S'-wave quarkonia Importance of the gluon fragmentation mechanism has been 
tested also in inclusive J/ip production in Z° decay Q]. There are still two open problems in 
the field. One is the polarization of prompt Jjjb at the Tevatron The other is the cross 

section for exclusive e + e~ — > J/ip + rj c [lllQ and e + e~ — > J/^ + cc + X at S-factories [3, 
Q]. However, the color-octet mechanism in NRQCD has been tested successfully in various 
ways. Comprehensive reviews on NRQCD phenomenology can be found in Ref. |15]. In 
the rest of this paper, we shall restrict our discussion to the gluon fragmentation into a 
color-octet spin-triplet QQ pair evolving into a S-wave heavy quarkonium. 

In the NRQCD factorization formalism the fragmentation function D(z,fi) for a par- 
ton splitting a heavy quarkonium is expressed as a linear combination of NRQCD matrix 
elements, which can be regarded as phenomenological parameters. Here, z is the momentum 
fraction of the final hadron relative to the decaying parton and \x is the hard-scattering scale 
of the process. Corresponding short- distance factors depend on z and are calculable in per- 
turbation theory. Most of the phenomenologically relevant short- distance factors have been 
calculated to leading order in a s . They all begin at order a 2 s or higher, with the exception of 
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the color-octet 3 Si term in the gluon fragmentation function, which begins at order a t 



a. 



The color-singlet 3 Si channel is suppressed because the short- distance factor begins at order 



a t 0. Since the color-octet 3 Si term dominates, the high-p^ gluon fragmentation phenom- 
ena in heavy quarkonium production, the next-to-leading order correction of order a 2 s to 
this term is particularly important. Unfortunately, two available results for the color-octet 



i S\ term disagree with each other 



Therefore, it is worth while to calculate this 



important function in an independent way. Since both previous calculations employed the 
Feynman gauge, we shall present our results in the light-cone gauge. 

The light-cone gauge is a physical gauge where the gluon field has vanishing light- 
cone projection A ■ n = 0, where n is an arbitrary light-like(n 2 = 0) vector appearing in 
the gauge-fixing term in the QCD Lagrangian. Derivation of the Altarelli-Parisi evolution 
of parton densities 3] is one of the best examples of the use of the light-co ne g auge. In 
a light-cone gauge, in which the fragmentation function was originally defined [l9|, there is 
a great simplification in QCD calculations. The eikonal line as well as the ghost decouples 
from the gluon, since the coupling, proportional to n M , is orthogonal to the gluon propagator. 
However, one draw-back in higher-order calculations is existence of the spurious pole 1/k-n 
in the gluon propagator 

k^n u + n^k u 



k 2 + 



-g 



hi/ 



k ■ n 



0" 



(1) 



where k is the momentum of the gluon. One should be very careful in dealing with the pole 
in evaluating loop integrals. As conventional methods, Cauchy-Principal- Value (CPV) or 



Mandelstam-Leibbrandt (ML) prescriptions 2(J |2l| have been used for a long time. The 



two prescriptions follow from different canonical quantizations, light-front (LF) and equal- 

ri 

time (ET), respectively [22[. It turns out that LF quantization (CPV) is untenable beyond 
the tree level as it conflicts with causality: the resulting theory cannot be renormalized. 
The ML prescription is causal; renormalization has been proved not only for the Green 
functions at any perturbative order but also for composite operators j^j]. All previous 
calculations beyond-leading order had to emp loy some prescription for the spurious pole. 



The CPV prescription is used in Refs 
leading order |2jj and the next-to-leading order [2j 
prescriptions can be found in Refs. 3(J, [3jJ, |3_ 



The ML prescription is employed in the 
Comprehensive reviews on these 



In our light-cone gauge calculation, we introduce a new method in determining the spu- 
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rious pole appearing in our calculation for the fragmentation function. We first employ a 
light-cone gauge to express the fragmenation function in terms of one-loop scalar integrals 
using the gauge-invariant definition of Collins and Soper [3^. However, we do not impose 
the sign of ie in the spurious poles and keep the scalar integrals from evaluation unlike 
the conventional ways. Our light-cone-gauge result for the fragmenation function is same 
as that for the Feynman gauge, which is guaranteed by the gauge-invariant Collins-Soper 
definition. Therefore, the scalar integrals involving spurious poles can be determined by 
comparing with the Feynman-gauge result. The spurious poles are finally identified with 
poles coming from the eikonal-line contribution in the Feynman gauge. Because all the 
poles are well-defined in the Feynman gauge, the spurious-pole contributions in the light- 
cone gauge are completely determined based on gauge invariance. Our result agrees with a 
previous calculation of Braaten and Lea Q m the Feynman gauge, but disagrees with the 
other previous calculation 16J. 



II. COLLINS-SOPER DEFINITION AND LIGHT-CONE GAUGE 

The fragmentation function D 9 ^h{z, fi) gives the probability that a gluon produced in 
a hard-scattering process involving momentum transfer of order /1 decays into a hadron H 
carrying a fraction z of the gluon's longitudinal momentum. This function can be defined 
in terms of the matrix element of a bilocal operator involving two gluon field strengths in a 
light-cone gauge Q. InRef. 0, 

Collins and Soper introduced a gauge-invariant definition 
of the gluon fragmentation function that involves the matrix element of a nonlocal operator 
consisting of two gluon field strengths and eikonal operators. One advantage of this definition 
is that it avoids subtleties associated with products of singular distributions. The gauge- 
invariant definition is also advantageous for explicit perturbative calculations, because it 
allows the calculation of radiative corrections to be simplified by using the Feynman gauge. 

The gauge-invariant definition of Collins and Soper for the gluon fragmentation function 
for splitting into a hadron H is j^] 

' n UJV-2 r+oo 



D ^ H (z, a) = — — / dx e 

x (01^(0)^(0-)^ V H(zk+fi±) £(x-) ba G+»(Q + , x~, ± )|0) . (2) 

The operator S(x~) in (J2J) is an eikonal operator that involves a path-ordered exponential 



of gluon field operators along a light-like path: 



Six )ba = Pexp 



+ig / dz-A + (0 + ,z-,0 ± ) 



(3) 



- ba 

where A^(x) is the matrix-valued gluon field in the adjoint representation: [y4^(a;)] ac = 
if abc A%(x). The operator Vh( p +, p± ) in Eq. (J2J) is a projection onto states that, in the asymp- 
totic future, contain a hadron H with momentum p = (p + ,p~ = (m 2 H + p 2 ±)/p + ,p±), where 
mil is the mass of the hadron. The hard-scattering scale p in Eq. (J2J) can be identified with 
the renormalization scale of the nonlocal operator. The prefactor in the definition (J2|) has, 
therefore, been expressed as a function of the number of spatial dimensions N = 3 — 2e. 
This definition is particularly useful when we use dimensional regularization to regularize 
ultraviolet divergences. If the production process of the hadron H can be described by 
perturbation theory, one can use the definition jj2J) to calculate the fragmentation function 
D 9 ^h{z, p) as a power series in a s . In Ref. complete sets of Feynman rules for this 
perturbative expansion for quark and gluon fragmentation functions are given. Since Ma 
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first used the definition |3J], the method became popular in quarkonium physics 
because of various reasons listed above. By inserting the eikonal operator (jSJ), the operator 
consisting of two gluon fields with different locations becomes gauge invariant. At higher 
order in a s , there are numerous diagrams which have gluons coupled to the eikonal lines. 
In the light-cone gauge, the contribution of the eikonal operator disappears since the gluon 
decouples from the eikonal line. Therefore, there is a great reduction in the number of 
Feynman diagrams. On the other hand, the spurious pole contribution of the gluon propa- 
gator appears in the light-cone gauge. However, the gauge invariance of this definition (J2J) 
provides the gauge transformation of the eikonal line contribution in the Feynman gauge 
into the spurious pole contribution in the light-cone gauge. By comparing the final results 
for the gauge- invariant quantity D 9 ^h{z, p) from the two gauges, the spurious pole coming 
from the gluon propagator in the light-cone gauge can be fixed unambiguously. 



III. PERTURBATIVE CALCULATION 



In this section we perform the next-to-leading-order calculation for the fragmenation 
function for a gluon splitting into a color-octet spin-triplet QQ pair of order a 2 s in a light- 
cone gauge. Although we carry out the calculation using a different gauge, we use the same 
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FIG. 1: Leading order Feynman diagram for g — > QQ. 

conventions as those presented in Ref. We do not reproduce the description on the 

theoretical background of the fragmentation function for heavy quarkonium production in 
NRQCD factorization formalism which is well explained in Ref. [l7|. Based on the NRQCD 
factorization formalism [4], the fragmentation function is written in a factorized form |l7|: 



Da-.H(Z) 



[(N - l)d T {z) + d L {z)\ W( 3 5r)) 



(4) 



where d,T and d^ are the short- distance coefficients for the transverse and longitudinal con- 
tributions and (Og( 3 Si)) is the color-octet 3 Si matrix element defined in Ref. ^J. 

There is only one lowest-order diagram in both Feynman and light-cone gauge, which is 
shown in Fig. ^ The circles connected by the double pair of lines represent the nonlocal 
operator consisting of the gluon field strengths and the eikonal operators in the definition 
(J2J). The momentum k = (k + , k~, k±) flows into the circle on the left and out of the circle 
on the right. The cutting line represents the projection onto states which, in the asymptotic 
future, include a QQ pair with total momentum p = (zk + ,p 2 /(zk + ), 0j_). The appearance 
of the diagrams for both gauges is the same in this order, since the circle should emit a 
gluon. With the Feynman rules of Ref. [33] and following the method of extracting the 



short-distance coefficients of the fragmentation function in Ref. 
order-a s terms in the short-distance functions dr(z) and di(z) as 

,2c 



d 



(LO). 



z) 



7ra s p, 



d 



(LO), 



8N(N - l)m 3 Q 
0. 



6(1 -z) 



17] . we can read off the 

(5) 
(6) 



We have neglected the relative momentum of the heavy quark in the QQ rest frame so that 
the invariant mass of the pair is p 2 = 4m 2 Q. The LO results (jSJ and © agree with previous 







FIG. 2: The Feynman diagrams of order a 2 for g — > QQ with final states. There are additional 
contributions from the complex-conjugate diagrams. 



calculations in the Feynman gauge [17j,|37[. Note that there is no reason for worrying about 
the spurious pole in this leading-order calculation. 

The Feynman diagrams for the fragmentation function for g — > QQ at order a 2 s consist 
of virtual corrections, for which the final state is QQ, and real-gluon corrections, for which 
the final state is QQg. The diagrams with virtual-gluon corrections to the left of the cutting 
line are shown in Fig. |2l The black blob in Fig. Efa) includes the vertex corrections and 
propagator corrections shown in Fig. El In the Feynman gauge, only the diagram in Fig. Ef b) 
vanishes, because the gluon attached to the eikonal line gives a factor of On the other 
hand, all the diagrams except for Fig. E(a) vanish in the light-cone gauge. If we use the 

38| | . we can simplify the structure of the 



threshold-expansion method of Braaten and Chen 
expression without employing the projection method. With the threshold expansion, we 
can keep the full structure of color and spin. Here we utilize the dimensionally regularized 
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(d) 



FIG. 3: One loop correction diagrams for g* — > QQ. 



threshold expansion method of Braaten and Chen [37, 39]. With the Dirac equation and the 
usual methods for reducing tensor integrals into scalar integrals, we factorize each virtual 
correction diagram into the leading order diagram in Fig. Q times a multiplicative factor. 
In the light-cone gauge, the ghost decouples since its coupling to the gluon is orthogonal to 
the gluon propagator (JTJ), so the gluon propagator correction factor shown in Fig. 3(d) does 
not have ghost contribution. 

The virtual corrections contribute only to the transverse short-distance function dx(z) 
defined in Ref. |l7|: 



d 



(virtual) 



4 LO) (^) x 2 Re 



A + II + 5Zr 



A 



(7) 
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where A is the vertex correction factor. The expressions 5Zq and IT are given by 

.2e r 

(2 - N) I ad + P 2 Iadd + 2p ■ nI BCD 
(4 - N) I ad + P 2 Iadd 



= ,16^ 



, I67ra s /j, 



2c 



n LC = -z6vra s /i 2e 



I! = — i 6%a s /j, 



2c 



1 2nt ( 1 



d+ iV 3 V iV 



/ab - 8p ■ uIabc 



lAB, 



(8) 
(9) 
(10) 

(11) 



where the superscripts F and LC stand for Feynman gauge and light-cone gauge, respectively 
The scalar integrals appearing in Eqs.([S |) - (jll|) are of the form 

d N+1 l 1 



Iab- 



(2ir) N + 1 AB ■ ■ ■ ' 

where the denominator AB ■ ■ ■ can be a product of 1, 2, 3, or 4 of the following factors: 



(12) 



A = I + ie, 

B = (I - p) 2 + ie = I 2 - 21 ■ p + Atuq + ie, 
C = (p — I) ■ n + ie, 

D = (I- p/2) 2 -m 2 Q + ie = l 2 -l-p + ie. 
The values for the integrals are in the Appendix of Ref. Q| except for 



(13) 
(14) 
(15) 
(16) 



47T 

tAAD = {A*)\2ml) 

i (An 
ADD ~ (4vr)2(2m2) I ^ 



r(i + e) 

em(l + 2e) 

£(l + e)_ 

em 



(17) 
(18) 



Note that for the light-cone gauge integrals Ibcd m Eq. (jHJ) and Iabc m Eq. (fTU|) have 
spurious poles l/(p — I) ■ n that require prescriptions. Those poles for the Feynman gauge 
are unambiguously defined as in Eq. (JT5j) . Therefore, we postpone evaluating the scalar 
integrals for the time being. If we use the ML prescription, each light-cone-dependent 



integral has ultraviolet (UV 
in the Appendix of Ref. 



and IR structures which are different from the values shown 
17| . Effectively, the ML prescription transforms a double pole 
into an IR pole and makes the integral satisfy naive power counting rules. Note also that 
Zq and IT are gauge dependent. The values in the Feynman gauge agree with well-known 
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ones that can be found, for example, in Ref. The result using the ML prescription 

is known only for the UV poles. We have full agreement in the UV poles if we use the 
ML prescription : the gluon propagator correction term is proportional to the QCD beta 
function as II = (33 — 2nf)a s /027re uv ) HQ and 5Z^ C = a s /(37re uv ) Q. All of them 
are listed, for example, in Ref. [42]. 

The contribution from the remaining diagrams shown in Fig. 2 (b)-(e), which have gluon 
couplings to the eikonal lines, is expressed as A. Their values are expressed in terms of 
one-loop scalar integrals: 



A 



LC 



2c 



9 7 



1 

N 



Iab^ + — -67) Iad-v Iaad 



+2p-n {9I ACD + Ibcd ~ 36/ A bc) 



2ira s fi 



2c 



1 



9 ( 1 + - ] I AB 



N + ^ - 47 ) I ad - P 2 Iaad 



A F = 
A LC = 0, 

A F = i 12 7ra s /i 2e I I A b - II AD + p-n (I A cd + Ibcd) 



(19) 

(20) 
(21) 
(22) 



The explicit value of the vertex correction factor A F in the Feynman gauge shown in Eq. (f2T)j) 
agrees with the result in Ref. H|. The UV dependence of the vertex correction factor A LC 
in the li ght- cone gauge shown in Eq. (|T9*|) agrees with the result using the ML prescription 



in Refs. 



41 



42| where only the UV contribution is given : Ak c = —5Zq C = — a s /(37reuv)- 



The integral Iaad has a Coulomb singularity as well as a logarithmic IR divergence due 
to the exchange of a gluon between the on-shell heavy quark and anti-quark. Dimensional 
regularization puts power infrared divergence like the Coulomb singularity to zero, so only 
the logarithmic IR divergence remains in the integral Iaad- Then the integral is effectively 
expressed by Iadd via the equation Iadd = (N — 4)Iaad- It is important to notice that 
various correction factors in the Feynman and the light-cone gauge involve different com- 
binations of the same scalar integrals. Straight-forward sums for both gauges produce a 
common result 



d 



(virtual) , 



d^°\z)^ Ret , 



7N- — + 51)I AD + 6n f [ I 



n ]Iab 



+ 18p-n {I ACD + Ibcd) + p 2 (81 add - Iaad)] }• (23) 
Thus the non-vanishing contributions from the gluon coupling to the eikonal line in the 
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Feynman gauge, A F , is simply distributed to other correction factors in the light-cone gauge 
via additional gluon propagator terms. 

Since gauge invariance holds for both the virtual and the real-gluon corrections sepa- 
rately, the equality of the virtual corrections in the Feynman and the light-cone gauge is 
a consequence of gauge invariance. As we commented previously, the light-cone dependent 
integrals in the Feynman gauge result have no ambiguities form spurious poles. On the 
other hand, we have not fixed the sign of the ie in the spurious pole of the integrals which 
are obtained in the light-cone gauge. Since we have found exact agreement between the 
two results in the two gauges, we may simply use the values obtained from the Feynman 
gauge calculation. Note that the integral Iabc disappears in Eq. ()23j). so the only light-cone 
dependent integrals that survive are Iacd and Ibcd- The values for these integrals are in- 
dependent of the sign of the ie in the definition of C in Eq. (|15p. The expansion of Eq. 
in e reproduces the result of Braaten and Lee 

L f 

,2 



4 virtual) (^) = 4 LO) (z)^'^ 



7i \ rrin 



3(1 - Q r ( l + e) + £(1 + + 177 - IQn, _ «_ + ^ + gh , g 



(24) 



2 euv^iR euv 18 2 

where (3q = (33 — 2n/)/6. 

The Feynman diagrams for the real-gluon corrections to the fragmentation function for 
g — ► QQ can also be calculated in both gauges. We draw the 5 left-half diagrams only, which 
must be multiplied by their complex conjugates to give a total of 25 diagrams. The real-gluon 
correction is a tree-level calculation. Therefore, there is no spurious pole problem. In the 
Feynman gauge, all 25 diagrams contribute, while only 9 diagrams in the light-cone gauge. 
In the latter gauge, diagrams^a) and^b) vanish. The real-gluon correction contribution is 
also gauge invariant. Employing either gauge, we reproduce the real-correction contribution 
given in Ref [u| before the phase-space integral is performed: 

d^(z) ~ ^ x 3 ° s (^X(l- 1 V T dx l — (25) 
{Z) - 8N(N - l)m% Trr(l-e) \m 2 Q J \ z(l-z)J J^ z)/z ' (25) 

where t — (1 — z)(zx + z — 1), x = 2q ■ p/p 2 , q is the final-state gluon momentum, and p is 
the QQ momentum. The final results for the real-gluon correction contribution of Braaten 
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FIG. 4: The Feynman diagrams of order a 2 for g — > QQ with QQg final states. There are a total 
of 25 diagrams, but only the left halves of the diagrams are shown. 



and Lee [17| are straight-forwardly reproduced: 



.(real), v 7Ta s /l 2e (X s ( 7T/i 2 



IN{N - l)m 3 Q 



rfi + e) 



7r \ m , 



x 



2euv e m 
6 /ln(l - z 

z 



z 



7TQ!, 



X 



1 - 2 

7T £ 



euv \il-zj+ 
+ 6(2 - z + z 2 ) ln(l - 2) 



1 — 2! 

+ + z(l-z] 



(27) 
(28) 



The infrared divergence cancels after summing the real and virtual correction contribu- 
tions shown in Eqs. (|24l) and ()27|) . Employing the MS scheme, a s and the operator are 
renormalized as in Ref. |l2|. After renormalization, the final answers for dr(z) and cIl(z) of 
Braaten and Lee [l^ are reproduced: 



d T (z,fi) 



7ra g (/i) 
48m^ 



5(1 + 



a s (//) 



7T 



A(n)8(l-z)+[\n^--l)P gg (z) 



-6(2- z + z 2 )\n(l - z) 



6 /On(l -z) 
1 - z 



(29) 
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where the coefficient is 



and Pggiy) is the gluon splitting function: 

z 1 — z 



7T" 



+ 81n2 + 61n 2 2 



6 



+ 



+ z(l-z) + ^- 5(1 - z) 
6 



(30) 



(31) 



16]. 



l-z) + z 

The transverse term dxiz) in Eq. (|2*^j) still disagrees with that of Ref. 1161 . Our final answer 
for the longitudinal fragmentation function is obtained by setting e — ► in Eq. (j2f 

o£(/x) l-z 



d L (z,fi) 



(32) 



The longitudinal term, dr (z) agrees with that of Braaten and Lee [ttJ as well as that of 
Beneke and Rothstein |43j. The dependence on the spectroscopic state of the produced 
quarkonium of this fragmentation function can be found in Ref. jl7j ]. 



IV. DISCUSSION 

We have calculated the next-to-leading order correction to the color-octet 3 Si term in the 
gluon fragmentation function in a light-cone gauge without using conventional prescriptions 
for the spurious pole. The gauge-invariant definition of the fragmentation function of Collins 
and Soper allows us to fix the ambiguities from spurious poles in the light-cone gauge 
by comparing with the result obtained in the Feynman gauge. Our result agrees with 
that of Braaten and Lee jl^l which disagrees with that of Ref. [la ]. The light-cone gauge 
considerably simplifies the calculation procedure for both the real and the virtual corrections. 
At least at one-loop level, the spurious pole problem can be resolved. This problem does 
not appear in the real corrections because they come from tree-level diagrams, but it does 
appear in the virtual corrections. We reduced the virtual correction in the color-octet 3 Si 
fragmentation function in the light-cone gauge to a linear combination of scalar integrals. 
After naive cancellations among the scalar integrals, ignoring the ambiguity from spurious 
poles, the correction reduces to scalar integrals that are independent of the sign of ie in 
the denominator k ■ n + ie. As a byproduct, the renormalization constants in the light-cone 
gauge were obtained at one-loop level. Their UV dependencies agree with the previous 
calculations within the ML prescription. They might be useful for other calculations, such 
as the next-to-leading order corrections to other fragmentation functions. 
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